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ABSTRACT 

We investigate the deflections of UHE protons by Galactic magnetic field(GMF) using four conventional GMF 
models in order to discuss the positional correlation between the arrival distribution of UHECRs and their sources. 
UHE protons coming from the direction around the Galactic center are highly deflected above 8° by the dipole 
magnetic field during their propagation in Galactic space. However, in bisymmetric spiral field models, there are 
dkections with the deflection angle below 1°. One of these directions is toward Centaurus A, the nearest radio- 
loud active galactic nuclei that is one of possible candidates of UHECR sources. On the other hand, UHE protons 
arriving from the direction of the anti-Galactic center are less deflected, especially in bisymmetric spiral field 
models. Thus, the northern hemisphere, not including the Galactic center, is suitable for the studies of correlation 
with sources. The dependence on model parameters is also investigated. The deflection angles of UHE protons 
are dependent on the pitch angle of the spiral field. We also investigate distortion of the supergalactic plane by 
GMF. Since the distortion in the direction around Galactic center strongly depends on the GMF model, we can 
obtain information on GMF around Galactic center if Pierre Auger Observatory finds the significant positional 
correlation around the supergalactic plane. 

Subject headings: cosmic rays — methods: numerical — IGM: magnetic fields — galaxies: general — 
large-scale structure of the universe 



1. INTRODUCTION 

The origin of ultra-high-energy cosmic rays(UHECRs) is one 
of the most intriguing problems in astroparticle physics. In or- 
der to reveal their nature and sources, UHECR observatories 
with larger exposures have been constructed. Now, the total ex- 
posure of Pierre Auger Observatory(PAO) already exceeds that 
of Akeno Giant Air Shower Array(AGASA), which had been 
the largest detector before the PAO era. Expectations for elu- 
cidation of UHECR sources are raised since PAO can detect 
more than 100 events per year above 4 x lO'^eV, which is an 
energy threshold that AGASA found small-scale anisotropy of 
observed UHECR arrival distribution from its 57 events. 

The discussions toward UHECR astronomy have been re- 
cently began. There seems to be two standpoints for UHECR 
astronomy. One is study of direct correlation between UHECR 
sources and observed UHECR arrival directions. This can be 
expected if the Galactic magnetic field(GMF) and extragalactic 
magnetic field(EGMF) are weak to deflect UHECR trajectories 
weakly. However, the nature of the GMF and EGMF is poorly 
known observationally and theoretically. Recently, several sim- 
ulations of the large-scale structure formation with magnetic 
field has describe d local mag r ietic s tructures (Sigl et al. 2003; 
iDolag et al.ll2005h . ISid etalJ(l2003h claimed that UHECR as- 
tronomy may not be possible if their source distribution model, 
magnetic structure and observer position are confirmed by fu- 
ture observation since UHE protons with 10^"eV are even de- 
flecte d 20° or mo r e by their structured EGMF. On the other 
hand, iDolag et"aLl ( l2005l) showed that most of highest energy 
cosmic rays are very weakly deflected by EGMF since strong 
magnetic field (~ /j,G) is highly localized at clusters of galax- 
ies. Therefore, the deflection angles of UHECRs are con- 
troversial. The understanding of the GMF is known better 
than that of EGMF, but poor, as briefly reviewed in section 



121 Number density of UHECR sources has also a crucial 
role. If the source number density is large, the number of ob- 
served events to unveil UHECR sources or source distribution 
is large. Recently, PAO reports the positional correlation be- 
tween the arrival directions of highest en ergy events and nearby 
AGN dPierre Auger Collaborationl2007l) . AGASA also reported 
small-scale anisotropy of UHECR arrival distribution within its 
angular resolution! Takeda et al. 1999). These results are one 
of corroborating evidences for this approarch. The UHECR 
source number density is constrained at 10"^ Mpc""' using 
small-scale anisotropy observed by AGASA, but this constraint 
has als o large uncert ainty due to small number of observed 
events iBlasi & de Marco, 2001 iKachekiess & Senukozll2005l: 



iTakami & Satol 2007ah . However, this uncertainty is wel l re- 
duced by nearly future observation dTakami & Satoll2b07al) . In 
any case, these uncertainty prevents our understanding of the 
positional correlation study at present. 

The other is purely statistical approach. On the arrival distri- 
bution of highest energy cosmic rays in data combined with re- 
sults of several observatories before the PAO era, medium-scale 
20°) anisotropy was found (Kachelriess & Semikoz 2006h. 
This angular scale corresponds a typical scale of event cluster- 
ings. The cumulative two-point auto-correlation function of the 
observed events can well reproduce that of mock data calcu- 
lated by Monte-carlo simulation assuming that UHECR source 
distribution t r aces ne arby large-scale structure within z — 0.02 
dCuoco et al.l l2007ah . They claimed that the medium-scale 
anisotro py reflects th e local large-scale structure as UHECR 
sources. ICuoco et a l. (2007b) suggested that such a statistical 
method is a powerful tool to reveal UHECR source distribution 
even if the GMF is conside red. This medium-sca le anisotropy 
was also found in PAO data jMoUerach et al.ll2007l) . 

Our recent studies have been based on the former standpoint. 
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In our previous study ( iTakami & Satol2007bh . we discussed the 
possibility that future observations of UHECR arrival distribu- 
tion can unveil large-scale structure of UHECR source distri- 
bution. We calculated the positional correlation between simu- 
lated arrival distribution of UHE protons and their source dis- 
tribution which reproduces the local structures, taking into ac- 
count their propagation in a structured EGMF with plausible 
strengths. Our EGMF model also constructed under simple as- 
sumptions so that it reproduces the observed local structures. 
Taking into account uncertainties on the EGMF strength and 
the source number density, we investigate the positional corre- 
lation in several numbers of their parameters. Our studies us- 
ing source distribution and EGMF model based on astronomical 
observations just predict the recent PAO result. We concluded 
that five year observation by R\0 can unveil nearby UHECR 
source distribution at a few degree scale if local number den- 
sity of UHECR sources is ^ 10"^ Mpc~^. More observation is 
needed if the local number density is larger The best indicator 
of the correlation is UHE protons above lO'^^eV. 

In the previous study, we considered EGMF as a first step to- 
wards the direct positional correlation study between UHECR 
arrival directions and source positions. The GMF was ne- 
glected. Ho wever, the GMF affects t h e arrival directions 
of UHE CRs ('Alvare z-Muniz et al.1 12OO2I: iTakami et all l2006t 
iKachelri ess et al. 200f~ Thus, as a next step, we should dis- 
cuss the deflections of UHECRs by the GMF. 

In this study, we calculate trajectories of UHE protons in the 
Galactic space and investigate their deflections by the GMF. We 
adopt the backtracking method for calculating the trajectories. 
Trajectories of particles with the proton mass and the charge 
of -1 injected from the Earth can be regarded as those of ex- 
tragalactic protons. We focus on UHE protons above 10'^ **eV, 
which is best indicator of UHECR sources. We also investigate 
the distortion of extragalactic sky by the GMF and discuss its 
effect to the correlation study. 

Compositions of UHECRs is essential for the positional cor- 
relation study. Heavier components of UHECRs are more 
strongly deflected by magnetic fields and disturb the positional 
correlation. One of observables for the composition measure- 
ments is the depth of shower maximum, X^a^^, which can be 
measured by fluorescence detectors. Its average value, (Xmax), 
is sensitive to UHECR composition. High Resolution Fly's 
Eye(HiRes) reported that compositions of cosmic rays above 
lO'^e V are domin ated by protons as a result of (X^ax) measure- 
ment (lAbbasi et al.. 2005) . Recent result by PAO i s compatible 
to the HiRes result within systematic uncertainty (jUngeretal] 
I2007h . Another observable is muon density from extensive 
air shower of UHECR, observed by the ground array. Recent 
studies of the muon content indicate some fraction of heavie r 
components in h ighest energy cosmic ravs lEngel et al.l ( l2007h : 
iGlushkov et al.l (l2007h . However, the interpretation of these 
two observables is dependent on hadronic interaction models, 
which include uncertainties at UHE energy region. On the 
other hand, the correlation with active galactic nuclei found by 
PAO indicates light composition. For the positional correlation 
study, proton component is a powerful tool. Thus, despite such 
uncertainties on UHECR compositions, we consider only pro- 
tons in this study. 

This paper is organized as follows: in section |2] we explain 
GMF models used in this study. In section[3j we investigate the 
deflection of UHE protons with 10'^ **eV in four GMF models 
introduced in section|2l Dependence on parameters in the GMF 



models is discussed in section |4] In section |5] we summarize 
results of the deflection and discuss effects of the deflection to 
the correlation study. 

2. GALACTIC MAGNETIC FIELD 

In this section, we briefly summarize current knowledge on 
the GMF and explain GMF models used in this stud y. For more 
details, see review articles (IValleell2004t lHanll200l . 

Faraday rotation measures of Galactic pulsars and extra- 
galactic radio sources indicate that the GMF has two compo- 
nents: large-scale regular fields with strength of few fiG and a 
turbulent or random field with comparable strength. The regu- 
lar component includes the disk and the halo components. 

The regular component in the Galactic disk has a pattern re- 
sembling that of the matter in the Galaxy. This spiral compo- 
nent could be well described by axisymmetric(ASS) or bisym- 
metric(BSS) field model. The BSS model has several number 
of field reversals and the ASS model does not. The halo field 
can be classified by parity at disk crossing. The odd parity:A 
(the even parity:S) represents that the spiral components in the 
halo fields above and below the Galactic plane are anti-parallel 
(parallel). From combination of the disk field models and the 
halo field models, four models are pro posed. It has discussed 
which model is favored. Vallee (2005*) reported that a concen- 
tric ring model (like an ASS model) with a single field reversal 
and S -type parity are p referred (ASS-S model). On the other 
hand, Han et al.l (1200 6') proposed that the spiral component has 
a bisymmetric structure with reversals on the boun daries of the 
spiral arms. In the halo field, CHan et al.lll997l[T99 9) suggested 
that A-type parity is preferable. This corresponds to BSS-A 
model. Both reports are based on observational results. How- 
ever, discrimination between the models is complicated by the 
presence of smaller scale irregularities in the magnetic field, 
as well as uncertainties in the theoretical modelling. Thus, we 
adopt all four combinations and investigate the deflections of 
UHE protons in the Galactic space. 

In the models, the radial and azimuthal components can be 
given by 

=B(r||,6')sin/9, Bg = B(r\\,9)cos p. (1) 

The r| | and 9 are the Galactocentric distance and azimuthal an- 
gle around the Galactic center respectively (9 is defined as in- 
creasing clockwise), p is the pitch angle of the spiral field in the 
neighbourhood of the solar system. We set p = -10°. The fie ld 
strength at a point (r| \,9)m the Galactic plane jStanevll 19971) is 
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Here, /3 = (tanp)"' = -5.67 and ro = 10.55 kpc is the Galacto- 
centric distance of the location with maximum field strength at 
I = 0° which can be expressed as ro = (/?0 + (i)exp(--|tan/7), 
where = 8.5 kpc is the distance of the solar system from the 
Galactic center and d = -0.5 kpc is the distance to the nearest 
field reversal from the solar system. Negative d means that the 
nearest field reversal occurs in the direction to the Galactic cen- 
ter b{r\\) is a radial profile of strength of magnetic field. The 
radial profile is modeled by 
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where Bo = 4.4/iG, which corresponds to 1.5/iG in the neigh- 
bourhood of the solar system. In the region around the Galactic 
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Fig. 1 . — Arrival directions of anti-protons with 10 eV injected from the Earth isotropically at 40 kpc from the Galactic center in the galactic coordinate. These 
can be regarded as arrival directions of extragalactic protons before modification by the GMF. 



center(r| | < 4 kpc), the field is highly uncertain, and thus as- 
sumed to be constant and equal to its value at r| | =4 kpc. The 
spiral field assumes to be zero for r\\ > 20 kpc. 

For the spiral halo field, we adopt an exponentially decrease 
with two scale heights (Stanev 1997) 

n< m ^^P^-^^ :0<z<0.5kpc 

where the factor exp(-3/8) makes the field continuous on z- 
The parity is represented as 

, / B(r||,0,z) : S type parity 

Wlh'P, z;- ..) : A type parity. 

In the GMF models with the A-type parity, the direction of 
the spiral field is reversed below the Galactic plane. As a re- 
sult, the spiral field is discontinuous at the Galactic plane. It 
is unphysical, but one of roughly approximate models reflect- 
ing observational results to support the A-type parity. Thus, we 
adopt such models to investigate the deflection of UHE protons. 
Note that GMF models with A-type parity predict symmetric 
trajectories of UHE protons about the Galactic plane even if the 
dipole field, introduced just below, is included. 

Near the solar system, the vertical component of magnetic 
field with strength of 0.2-0.3m G, direct ing toward the north- 
ern Galactic pole, is observed (i Han & Q iao 1994). Near the 
Galactic center, many non-thermal gaseous filaments perpen- 
dicular to t he Ga lactic plane with tens of /iG to mG have been 
discovered (lHanl l2Q07). These vertical magnetic fields indicate 
another regular component. In the dynamo theory, a dipole field 
is predicted with A-type parity as so-called AO mode. However, 
in this study, we assume that the z-component of magnetic field 
is a dipole field as 

Bx = -3 /iG sin 9 cos 9 cos ip/r' 

By = -3 /iG sin 6 cos 9 sin ip/ r' (6) 

B, = iic (1- 3 cos- 6l)/r^ 



in all four models. Here, 9 and ip are the zenith angle and the az- 
imuthal angle in the spherical coordinate centered at the Galac- 
tic center, respectively, /xg ~ 184.2 /iGkpc-' is the magnetic 
moment of the Galactic dipole which is normalized at 0.3/iG in 
the vicinity of the solar system. 

3. DISTORTION OF UHE SKY 

In this study, we investigate the deflection of UHE protons 
by the GMF. We adopt the backtracking method for calcula- 
tion of UHE proton trajectories in the GMF. We inject protons 
with the charge of - 1 (called anti-protons below) from the Earth 
isotropically, and follow each trajectory until the proton reaches 
a sphere of the Galactocentric radius of 40 kpc. The trajec- 
tory is regarded as that of a proton coming from extragalactic 
space. Their energy loss process is neglected since the energy 
loss length is typically much shorter than the propagation path 
length in the Galactic space. 

Figure [T] shows the velocity directions of anti -protons with 
jQi9.8gY jjjjgcfgjj from the Earth isotropically at 40 kpc from 
the Galactic center The maps are calculated in the cases of the 
ASS-S model(upper-left), the ASS-A model{upper-right), the 
BSS-S model(/ower-/e/f), and the BSS-A model{lower-right). 
Each point is regarded as an arrival direction of an extragalactic 
proton before modification by the GMF, e.g. its source direc- 
tion. 

All GMF models predict holes, where there are little ar- 
rival cosmic rays, with the radius of ~ 30° in the direction of 
the Galactic center Anti-protons injected from the Earth are 
strongly scattered around the Galactic center due to the strong 
dipole field. Thus, protons coming from the direction of the 
Galactic center do not arrive from extragalactic sources with its 
direction (see also Fig|3]l. 

In the cases of GMF models with the A-type parity, extra- 
galactic protons coming from the directions around the Galac- 
tic plain cannot also reaches the Earth. The main causes are the 
spiral field near the solar system with A-type parity. The tra- 
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Fig. 2. — Plots of the injection directions of anti-protons witli 10"^eV injected from tlie Eailli isotropically with colors which represent the deflection angle 
which the anti-protons experience during their propagation in the GMF in the galactic coordinate. Each point can be regarded as the arrival direction of a proton at 
the Earth with the deflection angle during its propagation in the Galactic space. 



jectories of protons arriving at the Earth are mainly affected by 
the spiral field in the Galactic plane and the dipole field around 
the Galactic center since strengths of these fields are relatively 
strong. Thus, anti-protons injected from the earth are, at first, 
deflected by the disk field near the solar system. Anti-protons 
injected to positive latitude are deflected to lower latitude in the 
case of injection toward the directions of the Galactic center, 
and to higher latitude in the case of injection toward the di- 
rections of the anti-Galactic center. In the both ASS models, 
the anti-protons keep on being deflected to the same directions 
due to no field reversal. Moreover, the A-type parity predicts a 
symmetric trajectories about the Galactic plane. Thus, the anti- 
protons are driven out from the Galactic disk region. On the 
other hand, the BSS models has several field reversals. Once 
the anti-protons pass over the reversal point, they start to be de- 
flected to the opposite directions to those just before. Therefore, 
the anti-protons are driven out less strongly from the Galactic 
disk region in the BSS-A than in the ASS-A model. In GMF 
models with S-type parity, the apparent separation does not oc- 
cur since the spiral field turns to the same direction above and 
below the Galactic plane. 

Next, we investigate the deflection angles of protons during 
their propagation in the Galactic space. In figure |2] the injec- 
tion directions of anti-protons with 10'^ '^eV are plotted with the 
deflection angles (shown as colors). In other words, a proton 
observed at a direction in the map has experienced deflection 
whose strength is represented as a color of the point. Note that 
the deflection angle is defined as a separation angle between 
the injection direction of an anti-proton and the direction of its 
velocity at 40 kpc from the Galactic center 

In any model, anti-protons injected to the direction of the 
Galactic center are strongly deflected by the dipole field near 
the Galactic center above 8°. The red regions with radius of ^ 
30° exists even if we consider anti-protons with lO^^ ^eV. Con- 
sequently, protons with the arrival directions near the Galac- 
tic center cannot positionally correlate with their extragalactic 



sources within a few degree. If there is n o dipole field, the 
red reg ions become much smaller as shown in lKachelriess et akl 
(|2007|) . We also discuss the effect of the dipole field in section 

HI 

In the both ASS models, the red regions also spread around 
the Galactic plane because of reflecting no field reversal as 
mentioned above. The region with relatively small deflection 
angles reflects the direction of the spiral field near the solar sys- 
tem. 

On the other hand, in the BSS models, anti-protons in- 
jected to the directions around the Galactic plane are less de- 
flected than those in the ASS models, but the deflection an- 
gles are 3° -6°. If anti-protons with 10^'"'eV are consid- 
ered, the dark green regions (representing the deflections with 
3° -5°) change into light-blue with l°-3°. The deflection 
angles except for the directions of the Galactic center are ap- 
proximately proportional to protons' energy. An intriguing 
structure is blue regions at around i£,b) ~ (-60° ,20°) in the 
both BSS models, and ~ (-60°, -20°) in the BSS-A models. 
The protons coming from these directions are very weakly de- 
flected since the deflections by the dipole field and the spiral 
field are balanced. In the former direction, there is Centaurus 
A(Cen A), which is the nearest radio-loud active galactic nu- 
clei(AGNs). That galactic coordinate is {l,b) - (-50°, 20°). 
Radio-l oud AGNs are one of str ong candidates of UHECR 
sources! T orres & Anchordoqu i 2004). PAO finds the positional 
correlation of the arrival d irections of highest energy event s 
with the direction of Cen A dPierre Auger Collaborationll2007h . 
If Cen A is one of UHECR sources and those models reflect 
the real universe, spatial correlation between Cen A and the 
arrival directions of highest energy is expected. If the dipole 
field is weaker than that in this calculation, the blue regions are 
shifted to a little lower longitude since the dipole field deflects 
anti-protons injected to low latitude and high latitude to higher 
longitude( Yoshiguchi et al. 2004). 

Next, we investigate the deflection direction of extragalactic 
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Fig. 3. — The arrival distributions of UHE protons with 10'"*eV predicted from UHECR sources arranged along the longitude lines and the latitude lines. 
These can be regarded as the distortion of the longitude and latitude lines. The longitude lines considered in this figure are £ = —\20° (yellow), —60° (red), 
0° (black), 60° (blue), 120° (light-blue), and 180° (green) respectively. The latitude lines considered in this figure are b = —60° (red), b = —30° (green), b = 0° (black), 
b = 30° (blue), and = 60° (light-blue) respectively. 
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Fig. 4. — The arrival distribution of UHE protons with 10"^eV predicted from UHECR sources arranged along the supergalactic plane. The colors represent the 
deflection angles of protons during their propagation in the Galactic space. 



protons and the distortion of extragalactic sky. Understanding 
global structure of UHECR deflections is important on the spa- 
tial correlation study. 

Figure [3] shows the distortions of extragalactic longitude 
lines(/e/f) and latitude lines{right). In left figures, we plot 
the injection directions of anti-protons with 10''^^eV whose 
velocity directions at 40 kpc from the Galactic center are 
£ = -120° (yellow), -6Q°(red), Q°(black), 6Q°{blue), \2Q°{Ught- 
blue), and \%Q)° {green) respectively. These lines are regarded 
as the arrival directions of extragalactic protons with lO'^'^eV 
from sources distributed along the longitude lines. 

In any model, there are many points with different colors 
around the direction of the Galactic center. These are the ar- 
rival directions of protons which are strongly deflected by the 
dipole magnetic field around the Galactic center. As also shown 
in FiglT] cosmic rays coming from the direction of the Galac- 
tic center arrive at the Earth from different directions in extra- 
galactic space. These cosmic rays hardly have information on 
the directions of their sources. About 5% of the events from 
the longitude lines of 60° and 300° and less than 3% of the 
events from the other three lines arrive from the direction of the 
Galactic center. Longitude lines near the Galactic center are 
distorted. The black lines are distorted by the dipole field. On 
the other hand. The red lines and the blue lines are also dis- 
torted, but by the spiral field just inside the solar system. Just 
inside the first field reversal interior t o the s olar system , there 
is strong spiral field (see also Fig. 1 in lAlvarez-Muniz & StanevI 
( l2006i) ). The distortion of the red and blue lines depend on the 
direction of this spiral field. The other lines seems almost not 
to be distorted. However, protons are actually deflected along 
the latitude line. The difference of the deflection between the 
GMF models are also apparent in the right figures. 

The right figures show the distortions of five latitude lines 
wifli b = -6Q°{red), b = -3Q°{green), b = Q°{black), b = 
3Q°(blue), and b = 6Q° {light-blue) respectively. The same 
method is used to describe the plots. The latitude lines are 
highly distorted, especially around the Galactic center. In 



the ASS-A model, the black line disappears since protons in- 
jected from extragalactic sources along the Galactic plane can- 
not reach the Earth, as also shown in FigUI In the BSS-A 
model, the black line almost disappears because of the same 
reason in the ASS-A model, but there are some black points 
around the Galactic center. However, these are deflected by the 
dipole field above 10° during their propagation, as we can see in 
Fig|2l Anti-protons injected from the Earth to positive latitude 
are deflected to higher latitude by the spiral field near the solar 
system in the case of the injection toward the direction of the 
anti-Galactic center Thus, the latitude lines with positive lati- 
tude around the anti-Galactic center are distorted to the lower 
latitude. The deflections are stronger in the ASS models than in 
the BSS models due to no field reversal, as discussed above 
and see Fig|2l On the other hand, the latitude lines around 
the Galactic center are distorted to the higher latitude in the 
northern galactic hemisphere, in addition to the deflection by 
the dipole field. These tendency are found in the right figures. 
In the southern galactic hemisphere, the deflection directions 
are unchanged in GMF models with S-type parity. In models 
with A-type parity, symmetric pattern about the Galactic plane 
is predicted. 

Finally, we investigate the distortion of the supergalac- 
tic plane. In recent years, PAO rejected that the hypothe- 
sis that the cosmic ray spectrum continues in the form of a 
power-law above an energy of lO'^'^eV with 6 sigma signif- 
icance dYamamoto et al.l I2007h . This fact, e.g. confirmation 
of Greisen-Zatsepin-Kuz ' m in(GZK) steepening ( GreisenI 19661 : 
IZatsepin & Kuz'rninlll966h . suggests astrophysical origins of 
UHECRs. Thus, the positional correlation between the arrival 
distribution of highest energy cosmic rays and the supergalactic 
plane is expected if near astrophysical objects are their sources. 
In the northern sky, the correlation between the arrival dis- 
tribution of cosmic rays above 4 x lO'^eV observed by four 
surface arrays(AGASA, Yakutsk, Haverah Park, and Volcano 
Ranch) and the supergalactic plane were pointed out with the 
chance probability from the uniform distribution of less than 
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Fig. 5. — Same as Fig|2] but in the cases of ASS-S and BSS-S models with an exponential radial profile of GMF strength(/e/?) and p = -8°(right). 



1% jUchihori et al.|l200Cl() . Recent PAO res ult indicates such 
correlation(|Pierre Auger C ollaborationl2007h . 

We investigate the distortion of the supergalactic plane by 
the GMF in order to revisit the correlation between the arrival 
directions and the supergalactic plane. Figure |4] shows the su- 
pergalactic plane{black points) and the arrival directions of pro- 
tons coming from the direction of the supergalactic plane(coZor 
plots). We plot the injection directions of anti-protons at the 
Earth(coZor points) whose velocity directions at 40 kpc from 
the Galactic center (black points) are the supergalactic plane. 
Each color shows the angle deflected by the GMF as shown in 
FigEl 

The supergalactic plane near the direction of the Galactic 
center is highly distorted by the spiral field inside the solar sys- 
tem and the dipole field around the Galactic center. Events from 
the direction of the Galactic center are highly deflected by the 
dipole field. In GMF models with S-type parity, the supergalac- 
tic plane is distorted to the opposite directions in the northern 
and southern galactic hemisphere. Thus, there is an intersec- 
tion between the supergalactic plane and distorted supergalactic 
plane. The direction of the intersection has smaller deflection 
angle. In GMF models with A-type parity, the supergalactic 
plane in the southern galactic hemisphere is distorted like the 
northern galactic hemisphere. In the BSS-A model, there are 
two intersection with the deflection angle of below 1 ° . 

The supergalactic plane near the direction of the anti- 
Galactic center is almost not distorted to the direction of lon- 
gitude, but distorted to the direction of latitude, as shown by a 
color of each point and also the right figures in Fig|3] In the 
ASS models, the deflection angles are larger due to no field re- 
versal. 

We summarize our results in this section. Protons coming 
from the direction of the Galactic center are highly deflected by 
the dipole field. Except for the direction of the Galactic center, 
the deflection angles in the BSS models are generally smaller 
than those in the ASS models because of the field reversals. Ex- 
tragalactic universe is distorted mainly along the latitude lines, 
which is reflecting local structure of the spiral field. In GMF 



models with A-type parity, extragalactic protons coming from 
the direction of the Galactic plane cannot reach the Earth. In 
the BSS models, protons coming from lower latitude are de- 
flected smaller than higher latitude since the latter protons do 
not experience strong modifications of their trajectories by the 
field reversals. In the BSS models, there are regions that the 
arrival directions are almost unchanged. One of the regions is 
near Cen A, which is the nearest radio-loud AGN which is one 
of strong candidates of UHECR sources. 

4. DEPENDENCE ON SEVERAL PARAMETERS 

In the previous section, we investigated the distortion of UHE 
sky by the GMF in four different GMF models with the same 
parameter set. However, the values of the parameters include 
some uncertainty. In this section, we discuss the dependence 
on such parameters. 

Figure |5] shows the same plots as Fig|2] using different pa- 
rameters in the ASS-S and the BSS-S model. A recent analy- 
sis of Faraday rotation measures suggests an exponential radial 
profile of magnetic field strength (Han et al, 2006). According 
to this analysis, we adopt the exponential profile supposed in 

iHan et aL IdSooi as 

b(r) = Boexp(^ ~^''~^^^ y (7) 

where r/, = 8.5 kpc is the scale radius. The two left figures in 
Fig|5]represents the maps of the deflection angles calculated in 
the GMF models with this exponential radial profile. The global 
structure is unchanged. 

The pitch angle also include some uncertainty. While 
we adopt p = -10° as the same value in prev ious studies 
(lAlvarez-Muniz et al.' '2002'; 'Tak ami et all |2006|). flie calcu- 
lations by Tinvakov & Tkachev ( 20021) and iProuza & Smidal 
(2003) adopted p = -8°. Here, we should check the dependence 
on the pitch angle. The two right figures in Fig|5]shows the de- 
flection angle maps calculated in the GMF models with the r"' 
profile with p = -8°. Bq is set to be 3.6/xG, which corresponds 
to 1 .5/xG near the solar system. 
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Fig. 6. — Same as Fig|3] but in the case of no dipole field. 



An interesting change is that the deflection angles in the di- 
rection of the anti-Galactic center are smaller than those in the 
models with p = -10° in the BSS-S model. The reason is the 
distance of a field reversal just outside the solar system. In 
the spiral field model we use, the nearest field reversal interior 
to the solar system exists at li = -0.5 kpc, which is a plausi- 
ble value based on many observations. On the other hand, the 
field reversals exterior to the solar system has weaker evidence. 
When p is changed, the field reversal points are changed except 
for the nearest one. In particular, the field reversal point just 
outside the solar system, which affects the deflection angles of 
protons arriving from the direction of the anti-Galactic center, 
become nearer than that in the case of p = -lQ°. Anti -protons 
injected to the anti-Galactic center are deflected to higher lati- 
tude by the spiral field in the neighbourhood of the solar system. 
When the anti-protons reaches the field reversal, the directions 
of the deflection are opposite and the trajectories are modified 
like their injection directions. If the field reversal is nearer, that 
modification is relatively strong. Thus, the deflection angles are 
smaller than those in the case of p = -10°. In the ASS-S model, 
similar tendency can be found. By definition, the spiral field 
is weak near points which correspond to the field reversals in 
the ESS model. The deflection angles of anti-protons injected 
to the direction of the anti-Galactic center are smaller since the 
region with weak magnetic field nears in the case of p = -8°. 

Finally, we investigate the effect of the dipole field. Figure 
|6] shows the same figure as Fig|3] but the dipole field is not in- 
cluded. As parameters of the spiral field, the same set as those 
in section [3] is adopted. Compared to Figl3] highly deflected 
events around the direction of the Galactic center disappear 
Anti-protons are less deflected since the spiral field is much 
weaker than the dipole field near the Galactic center However, 
because the magnetic strength of the spiral field reaches near 
10 /iG around the Galactic center, large deflections along the 
latitude lines are experienced by protons with their arrival di- 
rections of near the Galactic center. On the other hand, protons 
coming from the directions of the anti-Galactic center are al- 



most unchanged since the dipole field is very weak outside the 
Galactic center region. 

5. DISCUSSION & CONCLUSION 

In this study, we investigate the deflections of UHE extra- 
galactic protons with lO'^^eV in the Galactic space by calcu- 
lating their propagation in the four different GMF models using 
the backtracking method in order to discuss the correlation be- 
tween the arrival directions of UHECRs and their source posi- 
tions. 

Protons arriving from the direction of the Galactic center are 
highly deflected by the dipole field. Source positions of such 
protons are quite different from their arrival directions. Thus, 
the Galactic center region is not suitable for the correlation 
study. If the correlation is found in the direction of the Galactic 
center, we can obtain strong constraint on the GMF around the 
Galactic center. However, in the BSS models, there are direc- 
tions around the Galactic center that the arrival directions are 
almost unchanged. One of the directions is toward near Cen 
A, which is the nearest radio-loud AGN which is one of strong 
candidates of UHECR sources. If Cen A is confirmed as one of 
UHECR sources, the BSS models with dipole field and one of 
UHECR sources are confirmed. 

On the other hand, the deflection angles of protons com- 
ing from the direction of the anti-Galactic center are relatively 
small for the BSS models. Those become smaller by adopting 
smaller pitch angle, p = -8° since the field reversal exterior to 
the solar system is nearer. Thus, the anti-Galactic center region 
with 1^1 < 30° is best region for the correlation study. This re- 
gion is in the northern hemisphere, where PAO cannot detect 
UHECRs. 

We also investigate the distortion of the supergalactic plane 
by the GMF. The supergalactic plane near the Galactic center 
is highly distorted and the patte rn of the distortion is depen- 
dent on the spiral field model. Ilde et al.l (1200 ll) showed that 
we can check whether UHECR source distribution is correlated 
with the supergalactic plane by O(IO^) event detection all the 
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sky This claim is based on UHECR propagation in extragalactic 
space. If the positional correlation between the arrival distribu- 
tion and the supergalactic plane are found in the future, we can 
obtain information on the global structure of the GMF. On the 
other hand, the supergalactic plane near the anti-Galactic cen- 
ter is less distorted than that near the Galactic center although 
the arriving cosmic rays are deflected along the longitude line 
in the ASS models. In the northern hemisphere which con- 
tains the anti-Galactic center, the correlation between the arrival 
distribution and the supergalactic plane can be tested, which is 
weakly dependent on the spiral field model. 

As discussed above, the northern hemisphere is more suitable 
for the positional correlation study. The direction of the anti- 
Galactic center does not bother us with uncertainty of magnetic 
field around the Galactic center. In our previous study, we pre- 
dicted the positional correlation between UHECR arrival direc- 
tions and positions of local sources at a scale of 2° x 2° taking 
into account only EGMF( Takami & Satoll2007bl) . Adding the 
results in this study, we can predict the spatial correlation at a 
scale of - 4° for the BSS-S model with p = -8° and ~ 6° for 



the same with = -10° in the direction around the anti-Galactic 
center. If the ASS model is a real situation, regions where we 
can find the positional correlation in small angular scale are 
highly restrained. 

In the northern hemisphere. Telescope Array is un- 
der constructi on and will start observation in the near 
future! Fukushi maet al.ll20()7l) . PAO also project s its norther n 
site with larger exposure than its southern site(Nitz et al . 20071). 
Extreme Universe Space Observatory(JEM-EUSO), which can 
detect UHECRs from International S pace Station with ex - 
tremely large exposure, is also proiectedfebisuza ki et al.l2007l) . 
The southern site of PAO will strongly contribute to our un- 
derstanding on UHECR nature and sources. As next step for 
UHECR astronomy. These detectors, which can observe in the 
northern hemisphere, will be requested. 

The work of H.T. is supported by Grants-in-Aid for JSPS 
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Scientific Research provided by the Ministry of Education, Sci- 
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